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We calculate the finite-temperature shift of the critical wavevector Qc of the Pokrovsky-Talapov 
model using a renormalization-group analysis. Separating the Hamiltonian into a part that is 
renormalized and one that is not, we obtain the flow equations for the stiffness and an arbitrary 
potential. We then specialize to the case of a cosine potential, and compare our results to well- 
known results for the sine-Gordon model, to which our model reduces in the limit of vanishing driving 
wavevector Q = 0. Our results may be applied to describe the commensurate-incommensurate phase 
transition in several physical systems and allow for a more realistic comparison with experiments, 
which are always carried out at a flnite temperature. 
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I. INTRODUCTION 

The Pokrovsky-Talapo\ffl2] (PT) model describes a 
large variety of systems displaying a commensu- 
rate/incommensurate (C/IC) transition, ranging from 
vortex depinning in type-II superconductor^ to adsor- 
bate layers on crystal surfaceJ^ and q uantu m Hall bilayer 
systems under a tilted magnetic field .^^^^ This model is 
closely related to the sine-Gordon model,'S] with the ex- 
tra feature of a characteristic driving wavevector Q im- 
posed through the cosine term. In this work we present a 
functional renormalization group calculation of the finite- 
temperature corrections to the mean field results. 

The Hamiltonian for the PT model is 



HpT — 



1 2 

-ps(V0) -tcos{(j)- Qx) 



(1) 



where r = {x, y) and (Pr = dxdy. At mean-field level, 
one approximates the thermodynamic free energy by the 
Hamiltonian itself, evaluated at the field configuration 
that minimises the energy for given boundary conditions. 
This neglects the contributions of all fluctuations away 
from the minimum, therefore becoming less accurate as 
the temperature is raised and the entropic contribution 
of fluctuations increases. In this approximation, it is 
straightforward to compare the free energy of a config- 
uration with the field following the driving wavevector 
(commensurate phase) to that of one in which the field 
no longer follows the imposed Q. As Q increases, the 
presence of the stiffness term, suppressing deviations of 
the field from uniformity, makes it more and more costly 
to remain in the commensurate phase, until a critical 
Qc is reached at which the incommensurate phase be- 
comes energetically favourable. The aim of this article is 
to investigate the effects of thermal fluctutations on this 
critical Qc- 

We shall use a functional renormalization group (RG) 
scheme to study the model at finite temperatures. Our 
approach is as follows: We perform a simple transforma- 
tion which maps the PT model to a sine-Gordon model 



with additional terms depending only on the total topo- 
logical "charge" of the system and on the driving wavevec- 
tor Q. The RG transformation does not couple the sine- 
Gordon part to the Q-dependent part. Taking advan- 
tage of this, we renormalize the sine-Gordon part of the 
Hamiltonian and obtain a long-wavelength effective ac- 
tion, which we subsequently use to obtain the new value 
of the critical Qc- 

The main technical complication one faces in con- 
structing a systematic renormalization group transfor- 
mation for the sine-Gordon model is the inability to ex- 
pand the cosine term in Eq. (jlj in powers of the field 
and keep only a finite number of these. There are three 
main reasons for this complication: First, the periodicity 
of the cosine is crucial and is destroyed by any finite- 
order Taylor expansion. Second, we are interested in 
the two-dimensional follows from simple power 

counting,'^ polynomial interactions involving any power 
of the field are relevant in two dimensions. Finally, we are 
not interested in the fixed point of the RG transforma- 
tion, but in the actual values of the various parameters 
after integrating out the degrees of freedom that we are 
not interested in. 

Not expanding the cosine means that, in diagram- 
matic language, we ought to keep track of an infinite 
number of diagrams to one loop. This problem has al- 
ready been solved for the case of the Wilsonian RG by 
F. Wegner and A. Houghton in Ref. [10] where these au- 
thors derive the eponymous exact renormalization group 
equation. This equation is also the limit of an approxi- 
mate recursion relation first given by Wilson.'l^ Wilson's 
approximate recursion relation has been applied to the 
problem of critical wetting in Refs. |12I13] because the 
effective field theories used in studying critical wetting 
share with our problem its dimensionality, preventing the 
use of the more usual perturbation methods. In general, 
functional renormalization group approaches are useful 
for cases where there is an effective potential with a non- 
trivial functional dependence on the field, such as the 
cosine term in Eq. ([ij. 

The exact functional renormalization equation of Weg- 
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ner and Houghton rehes on a sharp, moving cutoff in 
momentum space. A sharp cutoff induces long-range in- 
teractions in real space and complicates the calculation 
of the flow of the stiffness ps in Eq. ([ij. This may 
be overcome by employing a smooth cutoff function j^'^ l ^^ l 
however, the resulting trajectories depend on the precise 
form of the cutoff. It has been showrP^l that the Wilso- 
nian RG approach (of which the Polchinksi RG is an 
example) suffers from strong scheme dependence even in 
the asymptotic regime. Various alternative formulations 
of RG transformations exist that do not suffer from this 
problem; one that has recently been applied to the sine- 
Gordon model is the functional renormalization of the 
effective average 

actioiJl3iS] (EAARG) in which a trans- 
formation is obtained, not for the Hamiltonian itself, but 
for the generating function for the 1-particle irreducible 
Green's functions. This is the RG scheme that we use in 
this article. 

The outline of this paper is as follows: In Sec. |ll] we 
describe the model and give a qualitative description of 



its behaviour. In Sec. Ill we explain the basic idea behind 
our approach before proceeding directly to the derivation 
(Sec. IV I and application to the PT model (Sec.|v} of the 
appropriate RG flow equations. Finally, we give a brief 
discussion of our results in Sec. IVII 



(/)(r) at some other point r' arbitrarily far from it either 
remains within a finite distance from 0o (smooth phase) 
or it does not but rather crosses over the maxima of V{(j)). 
Clearly, the system with t = is in the rough phasepB 
the usual RG analysis of the sine-Gordon model, Eq. ([iJ 
with Q — Q shows that, for given temperature T and 
Ps/T < I/Stt, there exists a below which the system is 
rough and above which it is smooth. For Ps/T > I/Stt it 
is always smooth. 

Consider now the case of finite Q. In the roughening 
picture, this corresponds to a potential V that depends 
on the position x; as x increases, the minima of the po- 
tential move to larger values of — 27rn -I- Qx. In other 
words, the potential is effectively "tilted". Thus, the po- 
tential part of the Hamiltonian tends to favour a cj) that 
increases with position and follows the potential, (p = (pn 
("commensurate phase"), while the gradient part favours 
a spatially constant cj) ("incommensurate phase"). The 
competition between ps and t leads to a transition be- 
tween the two states as, for example, t is varied. 

Notice that there are two separate effects here: one is 
the roughening transition (belonging to the Kosterlitz- 
Thouless universality class), which is already present 
when (5 = and the other is the commensurate- 
incommensurate transition, which appears only for finite 
Q. These two effects may be conveniently separated out 
as described in the next section. 



II. THE MODEL 



To qualitatively understand the features exhibited by a 
system described by the PT Hamiltonian given in Eq. 0, 
consider first the case Q = 0. It is then clear that, at 
mean field level, (j) will simply remain at one of the min- 
ima of the potential ¥{(/)) = —tcos{4>) at — 2nn with 
n = 0,±l,±2, ... 

On the other hand, consider the quantity 
(0(r)^(r')), = (0(r)^(r')) - (</>(r)) (</)(r')); for low 
enough temperatures, it is given by the 1-loop result 



('^(r)0(r')>, 



(2) 



in which ^ = y^ps/t, Kq{x) is a modified Bessel function 
and (. . .)^ indicates a cumulant. Since ^ ^ c» if i ^ 
and Ko{x) = ln(2/a;)-t-const-|-0(a;^ ln(a;)), the right hand 
side of Eq. (|2| diverges as i — > 0, i. e., there is no long- 
range order in the system. T his is merely an example of 
the Mermin- Wagner theorem.lSHSJ in addition, the Q = 
system exhibits a Kosterlitz-Thouless type transitioiJl^ 
on some line t{ps), which is again completely missed by 
a mean-field analysis. 

This transition is analo gous to the roughening transi- 
tion in interface problems .l^^^ In this analogy, the phase 
(j) corresponds to the height of an interface above a ref- 
erence level.SSI The phase in which ((/)(r)(/)(r'))^ remains 
finite as |r — r'| oo is called the "smooth" phase, 
while the one in which ((/)(r)0(r'))^ diverges is called the 
"rough" phase. Given the value of 0(0) = (po at some ar- 
bitrary point which we take to be the origin, the value of 



III. DETERMINATION OF Qc AT MEAN-FIELD 
LEVEL 



A. Separation of the Hamiltonian 

We begin by shifting to the new variable 6 — 
whereupon the Hamiltonian becomes 



with 



Hpt[0] — HsG + Hq, 



^Ps {Vef-t cos (9) 



and 



^jdxdJ = Y.<-^ 



' Qx, 
(3a) 

(3b) 
(3c) 

(3d) 



that is, the sum of the "charges" of all solitons present 
(per unit length); the integer i simply indexes the soliton. 

This form of the Hamiltonian is advantageous in that it 
consists of a sine-Gordon part, Hgc, which is independent 
of Q, plus the two terms in Hq, which do depend on Q. 
This last term is the essential difference from a simple 
sine-Gordon model. As we shall see in Sec. IV Hq is 
unaffected by an RG transformation. This will form the 
basis of our treatment of finite temperature effects. 



3 



B. Determination of critical Qc 

We shall next compute the critical Qc in the mean-field 
approximation; this will be straightforwardly extended to 
the renormalized case in Sec. |Vj 

To obtain the critical we notice that the transition 
from the commensurate to the incommensurate phase is 
signalled by the appearance of a finite soliton density 
with charge — We therefore calculate the excess en- 
ergy per unit area of a configuration with a single soliton 
(with charge -nP = —1) over that of one with no solitons, 
6* = 0; we denote this by Ai?. The part of the energy 
cost of a solitonic configuration due to Hsg is calculated 
in a standard way (sketched out in Appendix |A} to be 
-E-soi = 8y/ pst, hence AE = Eso\ — 2'iTpsQ. This quantity 
vanishes at the transition point, yielding a mean-field 
critical Qc oPI 



Qc 



TT V Ps 



(4) 



Note that Q^ diverges as 0, implying that the sys- 

tem remains in the commensurate state for all momenta; 
this is in agreement with the discussion at the end of 
Sec. |ll] according to which the creation of solitons (hence 
the transition to the incommensurate state) is caused by 
the stiffness overcoming the tendency of the phase field 
(j) to follow the minimum of the tilted potential. 



IV. EXACT RENORMALIZATION GROUP 
EQUATIONS 

Various schemes have been developed to study renor- 
malization group transformations of two dimensional 
field theories: In the theory of critical wetting, there 
have been several studies using this formulation of the 
RG, initially to first order in the potentiaP^ and then to 
all orders.li^ All these approaches have essentially used 
a local-potential approximation (LPA), in which the po- 
tential is allowed to change under coarse-graining, while 
the stiffness is not. The same method was later extended 
to a nonlocal model.l^Slln the same context there has also 
been work in which the gradient term is renormalized in 
an approximate way.E^ The LPA is generally thought to 
be applicable to wetting problems because the anomalous 
dimension is expected to be zero. 

In the present problem, on the other hand, it is 
clear that for vanishing Q one should obtain Kosterlitz- 
Thouless behaviour; furthermore, as discussed above, the 
quantity Q only couples to the number of solitons, which 
is conserved under the RG fiow. Thus, the appropri- 
ate RG must capture the Kosterlitz-Thouless type of be- 
haviour, which requires ps to fiow under the transforma- 
tion. 

As mentioned above, we shall use the effective average 
functional RG scheme introduced by WettericlJi^and ap- 



plied recently to the sine-Gordon model.li^For complete- 
ness, we outline the derivation of the exact fiow equation 
for the effective action before applying it to the sine- 
Gordon model; for more details, see Refs. I17I18I231 

From this point onwards, we will subsume the temper- 
ature into the parameters ps and t; that is, we use units 
in which the temperature kgT = 1. 

One begins by defining the bare action S[9] = H[d] 
and adding to it a piece 

A5h[0] = J^i?A(q)0(q)r(-q) 



where i?A(q) is called a regulator function (see be- 
low). One also adds a source term ASs[j,d] = 
S,(j(q)r(-q)-Hi*(-q)0(q)) and writes SA[j,e] = 
S[9] + ASnid] + ASs[j, 0]. The quantity 

WA[j]^log J V0expi-SA[j,O]) 

is then the generator of connected correlation function^ 
for the action S\. Its Legendre transform is 

f a[5] + Wa[j] = Eq fj*(-q)5(q) + j(q)^*(-q)) , 



where 9{q) = SW[j]/Sj*{—q) is the average of the field. 
We also define a new, related function 



TA[d]=fA[e]-ASR[0]. 



(5) 



Using well- knownP properties of W and F, and writing 
e = ln(Ao/A) where Aq is the initial value of the cutoff 
A, one finds 



deTA[0]\r- 



with T^^\0) indicating the n"* functional derivative of 
Fa. This is an exact result.E^It can be showrfl^that if 
RA{q) — > oo as A ^ cx) then Fa [6*] Sa[0]- fiuctuations 
about the mean-field solution are completely suppressed. 
Conversely, if i?A(g) ^ as A ^ then Ta[0] ^ T[e], 
so that the full generator of 1 particle irreducible (IPI) 
vertices is obtained. 

A full solution of Eq. ([6| for Fa would amount to com- 
puting all IPI functions of the system at some length 
scale ~ 1/A, including the full effects of fiuctuations. 
This is not a simple problem, and one must resort to 
approximations. We will take the form of Fa [6*] to be 



Tr[d,RA{q) Ff [0] + i?A(q) 



(6) 



d'r[-p,ie)iVey + Vie,e) 



A tedio us but straightforward computation leads 

^ j23l24l25l 
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(7a) 
(7b) 



t 




Figure 1: Flow diagrams for the RG Eqs. (10 I and (111 



where the operator V projects the function to its right 
onto the field-independent functional subspac^SU g^Yid 
-4(q) = {psQ^ + i?A(q) + V'-'^'^)-\ We take the cutoff 



function to be 



RAiq) 



(8) 



The parameter b controls the sharpness of the regulator 
function Ra in both wavevector and real space: for large 
b, i?A is local in wavevector space and long-range in real 
space; for b ^ 1, it is instead smooth in momentum space 
but sharp in real space. 

If V{e, 9) is restricted to its leading Fourier component 
V{e, 6) = —t{e) cos{d), the fiow equations are, to leading- 
order and after rescaling, 



det{e) - 2 



1 



AttPs 

iPsi^)) 



2-2/b 



n 



with n = bT{3 ~ 2/6)r(l + l/6)/(487r), reproducing 
the well-known leading- order fiow equations for the sine- 
Gordon model.!HEIl2lEa 

Including higher-order terms on the right hand side 
of Eqs. (7a I and (7b I is, in principle, straightforward; 



for example, expanding .A(q) in powers of V(t^ 0) and 
computing the integrals, one obtains 



F(2) (-l + 6)7r(f)'+^(\^(2))2 (_2 + 6)(-H-6)7r(^)"'" (F(2))3 



2bKp 



Ab^p^ sin(^) 



36V'^sin(^) 



(9) 



A similarly ungainly equation holds for d^psie). These 
expressions are greatly simplified if we make the choice 
& = 1 in the regulator function Eq. (|8|. For the poten- 
tial, the series on the right hand side of Eq. ^ may be 
summed (see Appendix IbI to yield 



A 

Airps 



In 1 



\/(2) 



Restricting the potential to the form V{£, 9) 



-t(e) cos(6') we obtain the fiow equation for i(e)ff 
1 



d^t{e) = 2t(e) - 



27rp,(e)i(e)/A2 



1- 1 1- 



A2 



(10) 

(see Appendi x p3| . In a similar way, and using the results 
of Appendix |B| to isolate the field-independent part, we 
obtain for the fiow of Ps(e) 



dePs{f) = 



247r(l -t2(e)/A4)3/2' 



(11) 



The fiow diagram corresponding to Eqs. (10 1 and (111 is 




pAe)l{kBT)Me)l{A^kBT) 



Figure 2: Evolution of A^(e) (dark full line) as a function of e. 
The dot-dashed line indicates the position of the minimum of 
^(e). The steep increase of ^(e) after its minimum is a result 
of our approximations and therefore unphysical. Notice that 
the minimum of A^(e) < 2-k. The initial values for this figure 
are Ps/keT = 0.06 and t/A'^ksT = 0.01. 



shown in Fig. ([T|. 

We now turn to the appUcation of the flow equations 
to the PT model. 



V. APPLICATION OF RG TO THE PT MODEL 
A. Calculation of Qc using the RG results 

To determine the scale at which we may stop inte- 
grating the flow equations and use mean field theory, 
one may use the scale-dependent correlation length (see 
Eq. (12)) ^(e) = ^/]^e)]t{7). Mean field theory applies if 

(e)<C 27r, while it is inapplicable otherwise. Thus, the 
appropriate e = £/ at which integration may be stopped 
may be located by integrating up to the point at which 
f (e) is a minimum; the mean-field approach of Sec. Ill 
may then be used to determine Q^. This is justified a 
posteriori if indeed A^(e/) < 2tt. This condition is satis- 
fied for all parameter values we have studied (see Fig. [5] 
for a representative example). 

How would the calculation described in Sec. IIIII be af- 
fected by the RG transformation? As we saw in Sec. |IV| 
Hq is invariant under the RG transformation, while Hsg 
is not, i. e., the parameters in HsG will change to Ps{() 
and t(e), respectively. To obtain the critical Qc, we 
notice that while Hq is unaffected by the RG, the en- 
ergy of a si ngle solit on now depends on e: we have 
Eso\{f) = 8-\/ps(e)t(e). This results in an energy differ- 
ence between the phase with no solitons and the phase 
with a single soliton given by AE — Eso\{f) — 2nps{Q)Q. 
Setting this equal to zero and solving for Q, as in sec- 
tion [ill] yields 




Figure 3: Demonstration of the asymptotic simple scaling 
of the trajectories, demonstrating their scheme independence 
(see text). The full line is t{e)/ ksT, the dashed line is 
Ps{e)/kBT and the dot-dashed line indicates the position of 
the minimum of ^(e) (see Fig. |2|. The steep increase of ^(e) 
after its minimum is a result of our approximations and there- 
fore unphysical (see text). The initial values for this figure are 
the same as for Fig. [2] ps/ksT = 0.06 and t/A^ksT = 0.01. 



The factor exp(— e) ensures that we are using physical 
(as opposed to rescaled) units. Eq. ( 12 I reduces to the 
correct mean-field expression, Eq. ([ij, for e = 0. 

In Fig. |2]we show a representative plot of the evolution 
of A^(e) with e. Evidently, the minimum is well below 
27r so that mean field theory is applicable to the renor- 
malized F. 

Let us now discuss the scheme-dependence of our cal- 
culation. In Ref. [16] it is shown that the trajectories re- 
sulting from the effective action functional RG scheme 
we use are scheme- independent, provided that the quan- 
tities t and ps flow as powers of the parameter e, i. e., 
Ps ~ exp((ipe) with some dp (and similarly for i). They 
term the region in which this occurs the "freezing region". 
Note that this scheme-independence does not hold in gen- 
eral for the case of Wilson-type renormalization. 

Fig. |3]is a log plot of the evolution of ps and t with e. 
The full line is t(e), the dashed line is Ps(e) and the dot- 
dashed line indicates the position of the minimum of ^(e) 
(see also Fig.[2|. Evidently, both Ps(e) and t{e) are in the 
freezing region at the values of e that we are interested 
in. This happens for all initial values of ps and t that we 
have checked. 



=(e) = exp(-e) 



7rp.(0) 



(12) 



B. Results and discussion 

In Figs. [4] (a) and (b) we show the ratio 7 = 
Qc(e)/Qc(0) determined by numerically integrating the 
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Figure 4: (Color online) The ratio 7 = Qc IQc'^ predicted 
by our RG analysis to its mean-field value. Figure (b) cor- 
responds to the lower left corner of Fig. (a). The strips on 
the right indicate meaning of shade. Notice the difference in 
scales between the two plots. White corresponds to 7 > 0.15. 



flow equations and applying the method discussed above; 
this corresponds to the ratio 7 — Q^'^ /Qc^^ ■ A darker 
colour indicates a larger decrease of the critical Qc due 
to thermal effects (see the insets and notice the differ- 
ent scales). The plot in Fig. [3] (b) is a zoomed- in part 
of Fig. [4] (a) (notice the range of the axes and also the 
changed color coding). 

For the purposes of this section, we will switch back 
to using natural units, by defining ps/ksT = ps and 
t/{kBTA^) = t, i. e., ps and t are the parameters in 
physical units (while t and pg are in units in which kgT = 
A = l). 

First let us fix p^ and T and vary i, i. e., fix pg and vary 
t. For increasing (decreasing) t, 7 increases (decreases), 
vanishing as t — 0. This occurs because, as described 
in Sec. [llj for vanishing t there is no long-range order 
even in the absence of Q. This effect is not present at 
mean-field level, hence 7 vanishes with decreasing t. 

Next, fix i and T and vary pg. From Fig. [i] 7 increases 
(decreases) for decreasing (increasing) pg or pig. 

Finally, consider fixing t and pg and varying the tem- 
perature T. This corresponds to fixed t/A'^pg or, in terms 
of Fig. |4] to moving on rays emanating from the origin 



with gradient t/pg. Increasing (decreasing) the temper- 
ature T corresponds to moving inwards (outwards) from 
the origin. As T is increased, a T will be reached at 
which the system will enter the rough phase at the ex- 
treme lower left corner of the diagram; thus, at high 
enough T, Qc vanishes for all pig and t. Equivalently, 
this may be viewed as a proliferation of solitons because 
the RG flow there is such that t{e) monotonically as 
e ^ 00, so that we also have £'soi(e) — > 0. This brings 
up the intriguing possibility of a purely temperature- 
driven commensurate-incommensurate transition in suit- 
able systems. For low enough temperatures, 7 — > 1 as 
expected. 



VI. CONCLUSIONS 

We have studied thermal effects on the commensurate- 
incommensurate transition point of the PT model using 
a renormalization-group approach. Our scheme relies on 
splitting the PT Hamiltonian into a sine-Gordon part, 
HgG, and a part depending only on the number of solitons 
present, Hq. We then derive a functional RG transfor- 
mation which acts on the sine-Gordon part while leaving 
the soliton part invariant. We are thus able to deter- 
mine the critical Qc at which the incommensurate phase 
eventually becomes stable, taking into account thermal 
effects. We find a general lowering of Qc compared to the 
mean-field result. Furthermore, there exists a regime in 
the pg-t plane in which Qc ^ even for finite t. This is 
due to the roughening transition, which is of Kosterlitz- 
Thouless type. 

Since real experiments are performed at finite temper- 
atures, the approach developed here may be applied to 
describe the C/IC transition in several physical systems. 
An application to the case of a Quantum Hall bilayer 
at total filling z/t = 1 in the presence of a tilted mag- 
netic field, motivated by recent experiments,'^^] will be 
described in a forthcoming publication.!^ 
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Appendix A: SOLITON ENERGY 

The excess energy of a single soliton may be calcu- 
lated as follows :_ We have E^oi = H[9] - H[0]; the 
solitonic profile 9 satisfies the extremization condition 
SHpT/66\g^g = 0. This condition gives an Euler- 
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Lagrange equation, the first integral of which yields 



dd 
2 V9a? 



= t (l — cos( 



(Al) 



Using this to eliminate the dependence of Hpt[6] on the 
cosine term and changing the variable of integration we 
obtain 



f'2TT 

Jo 



Finally, inserting Eq. (Al I into Eq. (A2| and integrating 
gives Eso\ = 8y/pst. 



Appendix B: SOME TRIGONOMETRIC 
IDENTITIES 



In this Appendix we display some of the calculations 



dx ' 



(A2) leading to Eqs. (|To| and ([iTl 

I 



It will be useful to know that 



cos™ (0) 



^((™")/2)cOs(0) + ..., 



odd 



' (w2) + 2^((™")/2)cos(20) + ... n 



even 



where the dots indicate higher-frequency Fourier components. 
We begin by expanding the integrand of Eq. (7a I for = 1, 



1 



1 r jPq 

2 J (27r)2 p,q2 + A2 + y(2) 



1 f (fq 



E 



(-T/(2))" 
2^ (2^)'f^i(Psg2+AT+' 



n \ A2 



1 (^^Y^ 



(Bl) 



(B2) 



Next, we find that 
In [1 + z cos (0)] ~ ^ — cos'' (0) 



= - (l - Vl^) cos ((/)) + 



(B3) 



where . . . indicates higher Fourier components. Eqs. (B2 1 
and (|B3| 



Finally, to restrict the flow of ps to the field- 
independent subspace, one needs the result 



2tt 



{l + zcos{(j))Y 2(i_^2) 



3/2 ■ 



(B4) 
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